An optical clock with neutral atoms confined in a shallow trap 
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We study the trap depth requirement for the realization of an optical clock using atoms confined 
in a lattice. We show that site-to-site tunnelling leads to a residual sensitivity to the atom dynamics 
hence requiring large depths (50 to 100 E r for Sr) to avoid any frequency shift or line broadening 



of the atomic transition at the 10 



level. Such large depths and the corresponding 



laser power may, however, lead to difficulties (e.g. higher order light shifts, two-photon ionization, 
technical difficulties) and therefore one would like to operate the clock in much shallower traps. 
To circumvent this problem we propose the use of an accelerated lattice. Acceleration lifts the 
degeneracy between adjacents potential wells which strongly inhibits tunnelling. We show that 
using the Earth's gravity, much shallower traps (down to 5 E T for Sr) can be used for the same 
accuracy goal. 

PACS numbers: 06.20.-f, 32.80.Qk, 42.50.Vk, 03.65.Xp 
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I. INTRODUCTION 

The control of the external degrees of freedom of 
atoms, ions and molecules and of the associated fre- 
quency shifts and line broadenings is a long standing issue 
of the fields of spectroscopy and atomic frequency stan- 
dards. They have been a strong motivation for the devel- 
opment of many widely spread techniques like the use of 
buffer gases Q , Ramsey spectroscopy 0, saturated spec- 
troscopy , two-photon spectroscopy^ , trapping and 
laser cooling 0,0, etc. 

In the case of ions, the problem is now essentially 
solved since they can be trapped in relatively low fields 
and cooled to the zero point of motion of such traps [j| . 
In this state, the ions are well within the Lamb-Dickc 
regime 1] and experience no recoil nor first order Doppler 
effect 5]. The fractional accuracy of today's best ion 
clocks lies in the range from 3 to lOx 10" 15 HEUEE!! 
with still room for improvement. The main drawback of 
trapped ion frequency standards is that only one to a few 
ions can contribute to the signal due to Coulomb repul- 
sion. This fundamentally limits the frequency stability of 
these systems and puts stringent constraints on the fre- 
quency noise of the oscillator which probes the ions ^3 ■ 

These constraints are relaxed when using a large num- 
ber of neutral atoms [13j for which, however, trapping re- 
quires much higher fields, leading to shifts of the atomic 
levels. This fact has for a long time prevented the use 
of trapped atoms for the realization of atomic clocks and 
today's most accurate standards use freely falling atoms. 
Microwave fountains now have an accuracy slightly be- 
low 10" 15 and are coming close to their foreseen ultimate 



limit which lies around 10~ 16 which is essentially not 
related to effects due to the atomic dynamics [lj, [l6| . In 
the optical domain, atomic motion is a problem and even 
with the use of ultra-cold atoms probed in a Ramsey- 
Borde interferometer > optical clocks with neutrals 
still suffer from the first order Doppler and recoil ef- 
fects 0, |2(1 . Their state-of-the-art accuracy is 
about 10" 14 [20f. 

The situation has recently changed with the proposal 
of the optical lattice clock [23] • The idea is to engineer 
a lattice of optical traps in such a way that the dipole 
potential is exactly identical for both states of the clock 
transition, independently of the dipole laser power and 
polarisation. This is achieved by tuning the trap laser 
to the so-called "magic wavelength" and by the choice 
of clock levels with zero electronic angular momentum. 
The original scheme was proposed for 87 Sr atoms using 
the strongly forbidden 1 Sq — 3 Po line at 698 nm as a clock 
transition |23j. In principle however, it also works for all 
atoms with a similar level structure like Mg, Ca, Yb, Hg, 
etc. including their bosonic isotopes if one uses multi- 
photon excitation of the clock transition [2J, |25j . 

In this paper we study the effect of the atom dynamics 
in the lattice on the clock performances. In ref. [22j, it is 
implicitly assumed that each microtrap can be treated 
separately as a quadratic potential in which case the 
situation is very similar to the trapped ion case and 
then fully understood |5| . With an accuracy goal in the 
10 -17 — 10~ 18 range in mind (corresponding to the mHz 
level in the optical domain), we shall see later on, that 
this is correct at very high trap depths only. The natu- 
ral energy unit for the trap dynamics is the recoil energy 
associated with the absorption or emission of a photon 



of the lattice laser, E r — 
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tor of the lattice laser and m a 
Sr and for the above accuracy goal the trap depth Uq 
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corresponding to the independent trap limit is typically 
Uq = 100 E r , which corresponds to a peak laser intensity 
of 25kW/cm 2 . 

For a number of reasons however, one would like to 
work with traps as shallow as possible. First, the resid- 
ual shift by the trapping light of the clock transition is 
smaller and smaller at a decreasing trap depth. The first 
order perturbation is intrinsically cancelled by tuning to 
the magic wavelength except for a small eventual tenso- 
rial effect which depends on the hyperfine structure of 
the atom under consideration. Higher order terms may 
be much more problematic depending on possible coin- 
cidences between two photon resonances and the magic 
wavelength j2^,|2^|. The associated shift scales as Uq |47j . 
The shifts would then be minimized by a reduction of 
Uq and its evaluation would be greatly improved if one 
can vary this parameter over a broader range. Second, 
for some of the possible candidate atoms, such as Hg for 
which the magic wavelength is about 340 nm, two-photon 
ionization can occur which may limit the achievable reso- 
nance width and lead to a frequency shift. Finally, tech- 
nical aspects like the required laser power at the magic 
wavelength can be greatly relaxed if one can use shallow 
traps. This can make the experiment feasible or not if the 
magic wavelength is in a region of the spectrum where no 
readily available high power laser exists, such as in the 
case of Hg. For this atom, a trap depth of 100 E r would 
necessitate a peak intensity of 500kW/cm 2 at 340 nm. 

When considering shallow traps, the independent trap 
limit no longer holds, and one cannot neglect tunnelling 
of the atoms from one site of the lattice to another. This 
leads to a derealization of the atoms and to a band struc- 
ture in their energy spectrum and associated dynamics. 
In section ITTTI we investigate the ultimate performance of 
the clock taking this effect into account. We show that 
depending on the initial state of the atoms in the lattice, 
one faces a broadening and/or a shift of the atomic tran- 
sition of the order of the width of the lowest energy band 
of the system. For Sr, this requires Uq of the order of 
100 E r to ensure a fractional accuracy better than 10 -17 . 

The deep reason for such a large required value of Uq 
is that site-to-site tunnelling is a resonant process in a 
lattice. We show in section IIVI that a much lower Uq 
can be used provided the tunnelling process is made non- 
resonant by lifting the degeneracy between adjacent sites. 
This can be done by adding a constant acceleration to the 
lattice, leading to the well-known Wannier-Stark ladder 
of states [27], . More specifically, we study the case 
where this acceleration is simply the Earth's gravity. The 
experimental realization of the scheme in this case is then 
extremely simple: the atoms have to be probed with a 
laser beam which propagates vertically. In this configu- 
ration, trap depths down to Uq ~ 5 E r can be sufficient 
for the above accuracy goal. 



II. CONFINED ATOMS COUPLED TO A 
LIGHT FIELD 

In this section we describe the theoretical frame used 
to investigate the residual effects of the motion of atoms 
in an external potential. The internal atomic structure 
is approximated by a two- level system \g) and |e) with 
energy difference ftw eg . The internal Hamiltonian is: 

Hi = huj eg \e)(e\. (1) 

We introduce the coupling between |e) and \g) by a laser 
of frequency to and wavevector k s propagating along the 
x direction: 

H s = Mlcos(u)t — k s x)\e)(g\ + h.c, (2) 

with Q the Rabi frequency. 

In the following we consider external potentials in- 
duced by trap lasers tuned at the magic wavelength 
and/or by gravity. The external potential H ext is then 
identical for both \g) and |e) with eigenstates |m) obeying 
H ex t\m) = hw m \m) (Note that |m) can be a continuous 
variable in which case the discrete sums in the following 
are replaced by integrals). If we restrict ourselves to ex- 
periments much shorter than the lifetime of state |e) (for 
87 Sr, the lifetime of the lowest z Pq state is 100 s) spon- 
taneous emission can be neglected and the evolution of 
the general atomic state 

\i>at) = £ al 4 \m, g)+a e m e"^+" m )* \ m> e ) (3) 

1 n 

is driven by 

*ftj^at> = (Hext + Hi + H.JIVat), (4) 

leading to the following set of coupled equations 

i< = ET elAm '' m * (m|e " lfcs£|m ' )a ™' (5) 

m' 
m' 

To derive eq. (jSJ) we have made the usual rotating wave 
approximation (assuming lj — Lu eg << uj eg ) and defined 

In the case of free atoms, H„ xt = 4-^- with hit the 
atomic momentum and m a the atomic mass. The eigen- 
states are then plane waves: \g, k) is coupled to |e, K + k s ) 
with A- - . r = oj — uj ea + hK/ - k ° -\- 7r -e-. One recovers the 
first order Doppler and recoil frequency shifts. 

Conversely in a tightly confining trap (m\e lksX \m' ^ 
to) << {m\e lksX \m) : and the spectrum of the system 
consists of a set of unshifted resonances corresponding to 
each state of the external hamiltonian. Motional effects 
then reduce to the line pulling of these resonances by 
small (detuned) sidebands . 
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III. PERIODIC POTENTIAL 
A. Eigenstates and coupling by the probe laser 

We now consider the case of atoms trapped in an op- 
tical lattice. As is clear from eq. (JSJ), only the motion of 
the atoms along the probe laser propagation axis plays a 
role in the problem and we restrict the analysis to 1D|48|. 
We assume that the lattice is formed by a standing wave 
leading to the following external hamiltonian: 



2m„ 



^(l-cos(2M)), (6) 
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FIG. 1: Band structure for two different lattice depth: Uo = 
2E r (left) and Uo — 10 E r (right). Each state |n, qo) is cou- 
pled to all the states \n , qo + k s ) by the probe laser. 



where hi is the wave vector of the trap laser. The eigen- 
states \n, q) and eigenenergies fku n ,q of H^ xt are derived 
from the Bloch theorem 29] . They are labelled by two 
quantum numbers: the band index n and the quasi- 
momentum q. Furthermore they are periodic functions of 
q with period 2ki and the usual convention is to restrict 
oneself to the first Brillouin zone q e] — kj.ki]. 

Following a procedure given in Ref. |28| a numerical 
solution to this eigenvalue problem can be easily found 
in the momentum representation. The atomic plane wave 
with wave vector k obeys 

^ t|K> = (l£ + t) |k> - t (|k+ 2kl) + |K - 2kl)) - 

For each value of q, the problem then reduces to the 
diagonalization of a real tridiagonal matrix giving the 
eigenenergies and eigenvectors as a linear superposition 
of plane waves: 



Hi xt \n,q) = Hu^ q \n,q) 



\n,q) 



E 



(8) 



with k 



l.q 



q + 2ifc;. For each value of q one obtains 
a discrete set of energies Hui^ q and coefficients C„ iKi<? , 
which are real and normalized such that J^i C% Ki q = 1 • 
In figures n and |21 are shown hu>^ q and Co )Ki , for various 
values of Uq. Except when explicitly stated, all numerical 
values throughout the paper are given for 87 Sr at a lattice 
laser wavelength 813 nm which corresponds to the magic 
wavelength reported in Ref. 30]. In frequency units E r 
then corresponds to 3.58 kHz. In figure|31 is shown the 



width (\^L q=kl 



a function of Uo in units of E r and in frequency units. 

Substituting (m\ — * (n,q\ and \m') — > \n',q') ineq. JSJ, 
the action of the probe laser is described by the coupled 
equations 



n,q 



l n,q+k s 



nn ,n 

-^— e " a n>,q+k 3 (9) 
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2 n',q 
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FIG. 2: C , Ki ,„ for two different lattice depth: U =2E r (left) 
and Uo = 10 E r (right). The bold vertical lines illustrate the 
case q = — ki/2. The dotted lines delimit the Brillouin zones. 
For a state \n = 0, q = aki) with a G] — 1, 1] the solid envelope 
gives the contribution of the plane waves Ift^afc; = a ki + 2iki). 



with n™> n 

uj eg + u^, iq - 0J n q+k 



nZiC n >,«,„C n , Ki , +ha and AJ." = u- 
n q+k e ■ As expected from the structure 
of the Bloch vectors in JSJ , the translation in momentum 
space e zkaX due to the probe laser leads to the coupling of 
a given state \n, q) to the whole set \n', q + k s ) (see figure 
^1 with a coupling strength f2™ *" and a shift with respect 
to the atomic resonance uu^, q+k — uu^ q . Both quantities 
depend on n, n' and q and to go further we have to 
make assumptions on the initial state of the atoms in the 
lattice. 
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FIG. 3: Lowest four band widths as a function of the lattice 
depth Uo in units of E r (left scale) and in frequency units 
(right scale). 
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FIG. 4: Left: Relative strength of the transitions to different 
bands (n = — * n) for an atom prepared in state \n — 0, q — 
—ki) (bold lines), \n = 0, q = —ki/2) and 171 = 0,5 = ki/2) 
(thin lines). Right: detuning of the first two sidebands for 
an atom prepared in state \n — 0,q — — fcj) (bold lines) and 
\n = 0, q = 0) (thin lines) in units of E r (left scale) and in 
frequency units (right scale). 
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FIG. 6: Expected resonances in the case where the first 
band is uniformely populated for Q = 10 Hz and Uo = 
20 E r , 30 E r , 40 E r , and 100 E r . The duration of the inter- 
action is such that the transition probability at resonance is 
maximized. 




FIG. 5: Shift of the "carrier" resonance in the first band for 
a lattice depth Uo = 10 E r . Left scale: in units of E r . Right 
scale: in frequency units. 



B. Discussion 

We first consider the case where the initial state is a 
pure \n, q) state. The strengths of the resonances ft™'" 
are shown in figure^Jfor the case n = and various values 
of q. At a growing lattice depth Q™'™ become indepen- 
dent of q and the strength of all " sidebands" (n' -n^0) 

asymptotically decreases as U ™'^ 4 for the benefit of 
the " carrier" [n' = n) . The frequency separation of the 
resonances rapidly increases with Uq (Fig. QJ. For Uq as 
low as 5E r , this separation is of the order of 10 kHz. For 
narrow resonances (which are required for an accurate 
clock) they can be treated separately and the effect of 
the sidebands on the carrier is negligible. If for example 
one reaches a carrier width of 10 Hz, the sideband pulling 
is of the order of 10~ 5 Hz. On the other hand, the "car- 
rier" frequency is shifted from the atomic frequency by 
u)^ q+kg — ujj l q due to the band structure. This shift is of 
the order of the width of the n th band (Fig. [5]and[3J|. 
It can be seen as a residual Doppler and recoil effect for 
atoms trapped in a lattice and is a consequence of the 
complete delocalisation of the eigenstates of the system 
over the lattice. The "carrier" shift is plotted in figure 
for the case n — and Uq = 10 E r . For this shift to be as 
small as 5 mHz over the whole lowest band, which corre- 
sponds in fractional units to 10~ 17 for Sr atoms probed 
on the 1 S'o — 3 Po transition, the lattice depth should be 
at least 90 E r (Fig. EJ. 

Another extreme situation is the case where one band 



is uniformly populated. In this case the "carrier" shift 
averaged over q cancels and one can hope to operate 
the clock at a much lower Uq than in the previous case. 
The problem is then the ultimate linewidth that can be 
achieved in the system, which is of the order of the width 
of the band and is reminiscent of Doppler broadening. 
This is illustrated in figure[|)]for which we have computed 
the expected "carrier" resonances in the case where the 
lowest band is uniformly populated, by numerically solv- 
ing equations JSJ. This was done for a Rabi frequency 
fi = 10 Hz and an interaction duration which is adjusted 
for each trapping depth so as to maximize the transition 
probability at zero detuning. We have checked that all 
resonances plotted in figure are not shifted to within 
the numerical accuracy (less than 10~ 5 Hz). However, 
at decreasing Uq the contrast of the resonance starts to 
drop for Uq < 40 E r and the resonance broadens progres- 
sively, becoming unusable for precise spectroscopy when 
the width of the energy band reaches the Rabi frequency. 
To get more physical insight into this phenomenon, let's 
consider the particular example of this uniform band pop- 
ulation where one well of the lattice is initially populated. 
This corresponds to a given relative phase of the Bloch 
states such that the interference of the Bloch vectors is 
destructive everywhere except in one well of the lattice. 
The time scale for the evolution of this relative phase is 
the inverse of the width of the populated energy band 
which then corresponds to the tunnelling time towards 
derealization (once the relative phases have evolved sig- 
nificantly, the destructive/constructive interferences of 
the initial state no longer hold) . The broadening and loss 
of contrast shown in figure can be seen as the Doppler 
effect associated with this tunnelling motion. 

The two cases discussed above (pure |n, q) state 
and uniform superposition of all states inside a band: 
J dq\n, q)) correspond to the two extremes one can ob- 
tain when populating only the bottom band. They illus- 
trate the dilemma one has to face: cither the resonance 
is affected by a frequency shift of the order of the width 
of the bottom band (pure state), or by a braoadening of 
the same order (superposition state), or by a combina- 
tion of both (general case) . In either case the solution is 
to increase the trap depth in order to decrease the energy 
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FIG. 7: External potential seen by the atoms in the case of 
a vertical lattice (Uo = 5E r ). An atom initially trapped in 
one well of the lattice will end-up in the continuum by tunnel 
effect. For Uo = 5E r the lifetime of the quasi-bound state of 
each well is about 10 10 s. 



width of the bottom band. 

In the experimental setup described in [3(J about 90 % 
of the atoms are in the lowest band and can be selected 
by an adequate sequence of laser pulses. The residual 
population of excited bands can then be made negligible 
(< 10~ 3 ). On the other hand, knowing and controlling 
with accuracy the population of the various \q) states in 
the ground band is a difficult task. The actual initial 
distribution of atomic states will lie somewhere between 
a pure state in the bottom band and a uniform super- 
position of all states in the bottom band. If we assume 
that the population of the \q) states in the ground band 
can be controlled so that the frequency shift averages to 
within one tenth of the band width, then a fractional ac- 
curacy goal of 10 -17 implies Uo = 70 E r or more. Note 
that due to the exponential dependence of the width of 
the ground band on Uo (see figure |3Jl the required lat- 
tice depth is largely insensitive to an improvement in the 
control of the initial state. If for example the averaging 
effect is improved down to 1 % the depth requirement 
drops from 70 E r to 50 E r . Consequently, operation of 
an optical lattice clock requires relatively deep wells and 
correspondingly high laser power, which, in turn, is likely 
to lead to other difficulties as described in the introduc- 
tion. 

Fortunately, the requirement of deep wells can be sig- 
nificantly relaxed by adding a constant acceleration to 
the lattice, as described in the next section. 



IV. PERIODIC POTENTIAL IN AN 
ACCELERATED FRAME 



Wannier-Stark states and coupling by the probe 
laser 



The shift and broadening encountered in the previous 
section are both due to site-to-site tunnelling and to the 
corresponding complete delocalization of the eigenstates 
of the lattice. As is well-known from solid-state physics, 
one way to localize the atoms is to add a linear com- 



ponent to the Hamiltonian [27J : adjacent wells are then 
shifted in energy, which strongly inhibits tunnelling. In 
this section we study the case where the lattice and probe 
laser are oriented vertically so that gravity plays the role 
of this linear component. The external hamiltonian is 
then: 



u 11 — 

n ext — 



2m a 



^(1 - cos(2k[£)) 



m a gx, 



(10) 



with g the acceleration of the Earth's gravity. This hamil- 
tonian supports no true bound states, as an atom initially 
confined in one well of the lattice will end up in the con- 
tinuum due to tunnelling under the influence of gravity 
(Fig. UJ. This effect is known as Landau-Zener tun- 
nelling and can be seen as non-adiabatic transitions be- 
tween bands induced by the linear potential in the Bloch 
representation |3l], H2, |33, |34( . The timescale for this ef- 
fect however increases exponentially with the depth of the 
lattice and for the cases considered here is orders of mag- 
nitude longer than the duration of the experiment • I n 
the case of Sr in an optical lattice, and for Uo as low as 
5 E r , the lifetime of the ground state of each well is about 
10 10 s! The coupling due to gravity between the ground 
and excited bands can therefore be neglected here. In 
the frame of this approximation the problem of finding 
the "eigenstates" of H^xt reduces to its diagonalization 
in a sub-space restricted to the ground band |35l |36( (we 
drop the band index in the following to keep notations as 
simple as possible). We are looking for solutions to the 
eigenvalue equation, of the form: 



H"\W n 



L ext\ vy m/ — 
\W m ) = 



(11) 



dqb m (q)\q) 



In eq. Ijlll) the \q) are the Bloch eigenstates of tl ext 
(c.f. section lirT)l for the bottom energy band (n = 0), m 
is a new quantum number, and the b m (q) are periodic: 
b m (q+2iki) = b m (q). After some algebra, eq. lfTT|) reduce 
to the differential equation 



h(w* - u™)b m (q) + im a gd q b m (q) = 



(12) 



where o A is the eigenvalue of the Bloch state \n = 0, q) of 
section lull Note that equations l|ll|l and 1)1 2|l only hold 
in the limit where Landau-Zener tunnelling between en- 
ergy bands is negligible. Otherwise, terms characterising 
the contribution of the other bands must be added and 
the description of the quasi-bound states is more com- 
plex [28|, [33 - In our case the periodicity of b m (q) and a 
normalization condition lead to a simple solution of the 
form 



,11 



,11 



mA, 



(13) 



b m (q) = 



-7,) 



with the definitions Wq 7 = 517 1-k dquj q , hA g 
m a gXi/2, and d q j q = lo 1 with 70 = 0. The \W„ 



FIG. 8: Wannier-Stark states in position (left) and momen- 
tum (right) representation for Uo = 5E r , Uo = 10 E r and 
Uo — 50 E r . Numerically we first compute the momentum 
representation (/t|Wm) = b m (K,)Co, K and then obtain the po- 
sition representation by Fourier tranformation. 



states are usually called Wannier-Stark states and their 
wave functions are plotted in figure |H1 for various trap 
depths. In the position representation \W m ) exhibits 
a main peak in the m th well of the lattice and small 
revivals in adjacent wells. These revivals decrease ex- 
ponentially at increasing lattice depth. At Uq — 10 E r 
the first revival is already a hundred times smaller than 
the main peak. Conversely, in the momentum repre- 
sentation, the distribution gets broader with increasing 
Uq. The phase shift between b m and fe m _i in 
b m (q) — e~™ q t kl b m -\{q), corresponds to a translational 
symmetry of the Wannier-Stark states in the position 
representation (x + Xi/2\W m ) = (a;|W m _i). The dis- 
crete quantum number m is the "well index" charac- 
terising the well containing the main peak of the wave 
function {x\W m ), and, as intuitively expected, the en- 
ergy separation between adjacent states is simply the 
change in gravitational potential between adjacent wells: 
HA g = m a g\i/2. 

Substituting (m\ — > (W m \ and \m') — > \W m ') in eq. 
(J3J) shows that the effect of the probe laser is to cou- 
ple the Wannier-Stark states to their neighbours by the 
translation in momentum space e lk " x , with the coupling 
strengths 



(W m \e ik °*\W m ,) = 




dn b* m (k + k s )b m > (K)C , K C , K +k. 



_(14) 

obtained from direct substitution of and (JBJ) [S0j |. 

Using the translational symmetry of the Wannier- 
Stark states it is easy to show that 

(W m \e ik ° x \W m ,) =e^ mk ^ kl (W \e lk ^\W m ^ m ). (15) 

From that property, equation l|14|) . and using b rn (n) = 
b* m (—K) (note that 7 g = 7- g ) one can then show that 

(W m \e lk ° x \W m+J ) = e^ k °/ k '(W m \e* k ° x \W m ^), (16) 

which is a useful result when studying the symmetry of 
coupling to neighbouring states (see next section). 

The differential equations JHJ), governing the evolution 
of the different states under coupling to the probe laser 



FIG. 9: Wannier-Stark ladder of states and coupling between 
states by the probe laser. 
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FIG. 10: Computed resonances when the initial state is a pure 
Wannier-Stark state. Left: Uo = 5E r , right: Uo = 10 E r . 
Both resonances are plotted for an effective Rabi frequency of 
the carrier = 10 Hz and for an interaction time of 50 ms. 



are then 



= ^T^^e-^e^--^ (17) 
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in which we have used and defined £7 m = 

n(W \e^ x \W m ) and A m = u - uj eg + mA g . 
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FIG. 11: Relative strength of the carrier |f2o/fi| 2 and of the 
first four sidebands |f2±i/f2| 2 and |f2±2/fi| 2 as a function of 
the lattice depth Uo- 
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FIG. 12: Frequency shift of the carrier as a function of the 
interaction duration in the case where the initial state of the 
atom is a coherent superposition of neighbouring Wannier- 
Stark states. Solid line: a%,(t = 0) = a g n+1 {t = 0) for all 
n. Dashed line: a 3 n (t = 0) = a 3 n+1 (t = 0)e iw/2 for all n. 
Dotted line: a B _ x {t = 0) = a 9 (t = 0) and a 9 n (t = 0) = for 
n —1,0. The shift is defined as the equilibrium point of 
a frequency servo loop using a square frequency modulation 
of optimal depth and computed for Uo = 5 E r and a carrier 
Rabi frequency Q.o/2n = 10 Hz. The interaction duration 
corresponding to a 7r pulse is tA g /2ir = 43.3. 



B. Discussion 

We now study the case where the initial state of the 
atom is a pure Wannier-Stark state. According to eq. 
(|17f) excitation by the laser will lead to a set of reso- 
nances separated by A g (see Fig. |S). In the case of Sr, 
A g /27T — 866 Hz and for the narrow resonances required 
for high performance clock operation, they are easily re- 
solved. The resonances obtained by first numerically in- 
tegrating l|14|) and then numerically solving (|17|) are plot- 
ted in figure ITUI for the cases Uq — 5 E r and Uo = 10 E r . 
They exhibit remarkable properties. First the "carrier" 
(which corresponds to the transition \ W m ) — > |Wm)) has 
a frequency which exactly matches the atomic frequency 
Lu eg . It also doesn't suffer from any broadening or con- 
trast limitation (provided the side bands are resolved) 
which would be due to the atomic dynamics. Second, the 
sidebands (|W m ) — > |W m ±»)) have a coupling strength 
which very rapidly decreases as Uq increases (see fig. I11H . 
In addition they are fully symmetric with respect to the 
carrier which results from eq. H16[) . and hence lead to no 
line pulling of the carrier. We have checked that the nu- 
merical calculations agree with this statement to within 
the accuracy of the calculations. This absence of shift 
and broadening remains true even for very shallow traps 
down to a depth of a few E r , the ultimate limitation 
being the lifetime of the Wannier-Stark states. This sit- 
uation is in striking contrast with the results of section 
IIIII in the absence of gravity. 

The system is more complex if the initial state of the 
atom is a coherent superposition of neighbouring wells. 
In this case off-resonant excitation of the sidebands will 
interfere with the carrier excitation with a relative phase 
which depends on the initial relative phase of neighbour- 
ing wells and on all the parameters of the atom-laser 
interaction (f2, to and the duration of the interaction). 
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FIG. 13: Frequency shift of the carrier as a function of the 
square modulation depth (see caption of Fig. 1121 . The cal- 
culation has been performed for Uo = 5E r , Q.q/2h = 10 Hz 
and an interaction time of tA g /2n = 43.5. The initial atomic 
state is the one corresponding to the dotted line in figure ED 



This interference leads to a modification of the carrier 
transition probability which is of the order of f2i/A 9 
(for the first, and most significant, sideband). For an 
interaction close to a n pulse, an order of magnitude of 
the corresponding carrier pulling is then f2if2o/A ff which 
can be significant. As an example for Uo = 10 E r and 
n /2n = 10Hz the shift is about 2 x 10~ 2 Hz, i.e. sev- 
eral times 10 -1 ' of the clock transition frequency. This 
shift is a priori all the more problematic as the initial 
atomic state is difficult to know and control accurately. 

We have numerically solved eq. (|17|l for various initial 
atomic states, lattice depths and interaction parameters 
to get a more quantitative insight of the effect. The re- 
sults are illustrated in figure ITSl for the case Uo = 5E r . 
A clear signature of the effect can be identified from its 
dependence on the interaction duration: the frequency 
shift oscillates with a frequency A g /2ir resulting from 
the A g term in A m _ m ' in p7jl. This provides a power- 
ful method for investigating site-to-site coherences in the 
lattice. More interestingly from a clock point of view, 
the shift becomes negligible for all interaction durations 
t such that t — (n + l/2)27r/A g . For these interaction 
durations the interference from the sidebands is symmet- 
ric for positive and negative detunings from resonance, 
leading to no overall shift. Since A g is extremely well 
known (potentially at the 10~ 9 level) this condition can 
be accurately met. Note that choosing such a value 
of the interaction duration does not significantly affect 
the contrast, as the two relevant timescales have differ- 
ent orders of magnitude in the narrow resonance limit 
(£! _1 >> A" 1 ), and therefore a range of values of t such 
that t = (n + l/2)27r/A s correspond to almost optimal 
contrast (e.g. all such values of t in figure IT^ . A more 
detailed study shows that the level of cancellation de- 
pends on the depth of the modulation used to determine 
the frequency shift (see caption of Fig. 1121) which results 
from a slight distortion of the carrier resonance. This 
effect is shown in figure ^3 which clearly indicates that 
the shift can be controlled to below 1 mHz even for a 
very shallow lattice depth down to Uo = 5E r . 
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V. DISCUSSION AND CONCLUSION 

We studied the trap depth requirement for the opera- 
tion of an optical lattice clock with a projected fractional 
accuracy in the 10~ 17 — 10 -18 range. We have shown that 
using a purely periodic potential necessitates a lattice of 
depth 50—100 E r limited by tunnelling between adjacent 
sites of the lattice. A possible way to vastly reduce this 
depth is to use gravity to lift the degeneracy between the 
potential wells which strongly inhibits tunnelling. Trap 
depths down to 5 — 10 E r are then sufficient to cancel 
the effects of the atom dynamics at the desired accuracy 
level. This will become even more important for future 
work aiming at even higher accuracies. Although very 
simple, gravity is not the only way to suppress tunnelling 
and other solutions, essentially consisting in a dynamic 
control of the lattice, are certainly possible |3Sl KM flof . 
They may prove useful if one wants to operate a lattice 
clock in space for instance. 

Throughout the paper we have not taken into account 
the dynamics of the atoms in the directions transverse 
to the probe beam propagation. Experimental imperfec- 
tions however (misalignemcnt, wavefront curvature, aber- 
rations) may lead to a residual sensitivity to this dynam- 
ics. If for example the probe beam is misaligned with 
respect to the vertical lattice by 100 /irad the transverse 



wavevector k± is about 10~ 4 k s and a modest transverse 
confinement should be sufficient to make its effect negligi- 
ble. Such a confinement can be provided by the gaussian 
transverse shape of the laser forming the lattice or by a 
3D lattice. The latter also leads to an interesting physi- 
cal problem depending on the relative orientation of the 
lattice with respect to gravity ^lj ■ 

Finally the well-defined energy separation between 
Wannier-Stark states and the possibility to drive tran- 
sitions between them on the red or blue sideband of 
the spectrum (section HVB() opens new possibilities for 
the realization of atom interferometers. This provides a 
method to generate a coherent superposition of distant 
states for the accurate measurement of the energy sepa- 
ration between these states. This can for instance lead to 
an alternative determination of g or h/m a |4lll43ll44l45| . 
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